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Abstract 

This paper deals with an optimal control problem related to a phase field system 
of Caginalp type with a dynamic boundary condition for the temperature. The 
control placed in the dynamic boundary condition acts on a part of the boundary. 
The analysis carried out in this paper proves the existence of an optimal control for a 
general class of potentials, possibly singular. The study includes potentials for which 
the derivatives may not exist, these being replaced by well-defined subdifferentials. 
Under some stronger assumptions on the structure parameters and on the potentials 
(namely for the regular and the logarithmic case having single-valued derivatives), 
the first order necessary optimality conditions are derived and expressed in terms 
of the boundary trace of the first adjoint variable. 
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1 Introduction 


In this paper we are concerned with an optimal control problem for a nonlinear phase 
held system of a standard form (cf. the monograph [6]), but with a possibly singular 
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double well potential, like in the logarithmic case (cf. the later fll.djl h and with a dynamic 
boundary condition for the temperature, in which also the time derivative of the boundary 
temperature plays a role and where the control variable appears in the external term 
(see fll.bjl h Let us now introduce and discuss the problem in precise terms. 

A rather general version of the phase held system of Caginalp type [7] reads as follows 

dt'd — A'd + X{(p)dt(p = 0 in Q := (0, T) x (1.1) 

dt^p — aA(p + W'((p) = 'dA(<p) in Q (1.2) 

where 12 is the domain where the evolution takes place, T is some hnal time, i? denotes 
the relative temperature around some critical value that is taken to be 0 without loss 
of generality, and is the order parameter. Moreover, A is a given real function, whose 
meaning is related to the latent heat, and a is a positive constant. Finally, W' represent 
the derivative of a double-well potential W. Typical examples for W are the regular 
potential 

AVreg(p) = r^(r — 1)^ , pGM (1.3) 

with two absolute minima located in 0 and 1, and the logarithmic potential 

'^iog{r) = ((1 + r) ln(l -|- r) 4- (1 — r) ln(l — r)) — ar^ , r G (—1,1) (1.4) 

where the coefficient a > 0 is large enough in order to kill convexity. The potential fll.3p 
is a shifted version of the usual classical potential given by r i—)■ |(r^ — 1)^ and precisely 
satishes our general assumptions given below, while fll.dp has a derivative which behaves 
singularly in the neighborhoods of —1 and 1. Generally speaking, the potential W could 
be just the sum W = j3 + where ^ is a convex function that is allowed to take the 
value -l-cxo, and ^ is a smooth perturbation (not necessarily concave). In such a case, f3 is 
supposed to be proper and lower semicontinuous so that its subdifferential {3 ■.= d[3 is 
well dehned and can replace the derivative which might not exist. Of course, W' has to 
be read as /3 -|- tt, where vr := vf', and equation fll.21) becomes a differential inclusion. 

Moreover, initial conditions like 19(0) = and (p(0) = <po and suitable boundary 
conditions must complement the above equations. As far as the latter are concerned, we 
take the homogeneous Neumann condition for (p , that is, 

dnP> = 0 on S := (0,T) x T (1.5) 

where T is the boundary of 12 and is the (say, outward) normal derivative. The 
position 01.51) is mostly the rule in the literature for the order parameter ip. Concerning 
the temperature i?, in order to address the boundary control problem, we choose the 
following dynamic boundary condition 

dn'& + rdt^Y + a('dr — iRw) = 0 on S (1.6) 

where := 'djp is the trace of d on the boundary and u is the control, which is supposed 
to vary in some control box Had- Moreover, in fll.bp . r is a positive time relaxation 
parameter, a is a positive constant and m is a nonnegative function dehned on T. Notice 
that, in fact, the control u can act just on the subset of T where m is positive. Thus, the 
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state system takes the following form 


dt'd — A-d + \{ip)dtip = 0 

in Q 

(1.7) 

dtip - a Alp + (3{lp) + 7r(<p) 3 '&\{(p) 

in Q 

(1.8) 

dn'& + rdt'dr + a ('dr — ’mu) = 0 and dn<p = 0 

on S 

(1.9) 

'd(O) = "do and <p(0) = <po 

on n. 

(1.10) 


The cost functional we consider depends on two nonnegative constants Ki and K 2 and two 
functions dq and on Q and hi, respectively. We want to minimize 

aw '’ol' + f ^ (I’ll) 

where m, 'd and (p vary under the constraint of the state system and u G Uad, where the 
control box Uad is dehned by 

Uad := {u e : Wmin <U< Umax a.c. ou S} (1.12) 

for some given bounded functions Wmin and Wmax- The analysis carried out in this paper 
shows the existence of an optimal control for a general class of potentials W = ^ + ^: 
indeed, for this purpose (3 is just assumed to be a general convex and lower semicontinuous 
function with minimum 0 attained at 0, that is, /3(0) = 0, which is physically reasonable. 
On the other hand, the derivation of the hrst order necessary optimality conditions can be 
made only in case of regular (like fll.dp i and singular (like 01.4^ 1 potentials. Linearized and 
adjoint problems are under our investigation and, subsequently, the optimality conditions 
can be expressed in terms of the adjoint variables (see the Theorem 12.91 stated in the next 
section). 

Let us mention here some related work. As far as we know, the contributions on 
optimal control problems for phase held models are quite a few and often restricted 
to the case of regular potentials, or dealing with approximations of the actual systems 
when the hrst order optimality conditions are discussed. In this respect, we point out 
the papers dH 12011211121] concerned with distributed control problems; we also refer to 
[211111I21II311IS112311201I2B] for diherent types of phase held models and other kinds of control 
problems. The main features of our paper are the study of a boundary control problem and 
the consideration of a dynamic boundary condition, in a very simple form: indeed, (ITTD 
is an affine condition involving the temperature and its time derivative on the boundary, 
with an external term carrying out the action of the control. About dynamic boundary 
conditions, also of nonlinear type and possibly involving the Laplace-Beltrami operator, 
let us quote the articles [HlISl ITIJlfmiT^ITTHn] . 

The paper is organized as follows. In the next section, we list our assumptions, state 
the problem in a precise form and present our results. The well-posedness of the state 
system, the regularity results and the existence of an optimal control will be shown in 
Sections |3] and lU respectively, while the rest of the paper is devoted to the derivation of 
hrst order necessary conditions for optimality, which are computed in the case of potentials 
that generalize fll.3p - fll.4p (some cases of more singular potentials being the subject of 
a future project of ours). The hnal result will be proved in Section El and it is prepared 
in Sections [5] with the study of the control-to-state mapping. Finally, the Appendix is 
devoted to the rigorous proof of an estimate that is derived just formally in Section |3l 
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2 Statement of the problem and results 

In this section, we describe the problem nnder stndy and present onr results. As in 
the Introduction, hi is the body where the evolution takes place. We assume hi C 
to be open, bounded, connected, and smooth, and we write |r2| for its Lebesgue measure. 
Moreover, F and dn still stand for the boundary of hi and the outward normal derivative. 


respectively. Given a finite final time T > 0, we set for convenience 

Qt := (0,f) X n and := (0,f) x F for every t G (0,T] (2.1) 

Q ■= Qt , and I] := ■ (2-2) 

Now, we specify the assumptions on the structure of our system. We assume that 

a, T E (0, +oo), m E L°°{Q,) and m > 0 a.e. in hi (2.3) 

/5 : M —)■ [0, +cxo] is convex, proper and l.s.c. with /3(0) = 0 (2.4) 

T , A : M —?■ M are functions and tt', X' are Lipschitz continuous (2.5) 

We set for convenience 

/3:=d/3, 7T :=Tf' and A := (2.6) 


and denote by D{P) and D{f3) the effective domains of (3 and /3, respectively. Moreover, 
(3°{r) is the element of (3{r) having minimum modulus for every r G D{f3) (see, e.g., [5|, 
p. 28]). It is well known that f3 is a. maximal monotone operator from M to M (see, 
e.g., O Ex. 2.3.4, p. 25]). Next, in order to simplify notations, we set 

V := H := ^^(fi), W ■.= {v E iF^(fl) : dnV = 0} 

as well as Hr := ^^(F) and W := i/^(F) (2.7) 

and endow these spaces with their natural norms. The symbol || • |]x stands for the norm 
in the generic Banach space A, while || • ||p is the usual norm in anyone of the H spaces 
on n, F, Q and S, for 1 < p < cxo, provided that no confusion can arise. Furthermore, it 
is understood that H is embedded in V*, the dual space of V, in the standard way, i.e., 
in order that {u, v) = uv for every u E H and v E V, where (•, •) denotes the duality 
product between V* and V. Finally, for v E L‘^{Q) the symbol 1 * n is the usual time 


convolution, i.e., 

{l*v){t):= I v{s)ds fortG[0,T]. (2.8) 

Jo 

At this point, we describe the state system. Given -do and ipo such that 

e 1^, <Po e IF and /3°((po) e H (2.9) 

we look for a quadruplet ('d,'dr,<F,0 satisfying 

E H\0,T;H)nL^{0,T;V) (2.10) 

idr E H^{0,T] Hr) and 'dr(^) ='d(^)|p for a.a. f G (0, T) (2-11) 

if E 1F^’°°(0, T; H) n 77^0, T; V) n L°°(0, T; W) (2.12) 

e G 7:°°(0,T;77) (2.13) 
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and solving the problem 


[ dt-dv + 

[ Vd 

■ Vn -7 

/ X{ip)dtipv + T 

/ dtiHr Vr + a (dp — mu] 

1 Up = 0 

Jn J 

'n 


m 

Jr Jr 


for every 

{v, Vr) G 17 X 

Vy such that up 

= n|p and a.e. in (0,T) 

(2.14) 

dtip — aAip 

+ e + 

7r{ip) = 

'dX{(p) and f > 

G (J{ip) a.e. in Q 

(2.15) 

dnip = 0 a 

.e. on 

S 



(2,16) 

m = ^0 

and 

<d(0) = 

ipo a.e. in 12. 


(2.17) 


Remark 2.1. The variational equation fl2.14jl is the weak formulation of equation fll.7j) 
and of the dynamic boundary condition contained in fll.91) . Let us notice that we can 
deduce both fll.71) and the hrst condition in fll.9p from 02.141) . Indeed, by writing 02.14p 
with an arbitrary v G Ro(fl) and ur = 0, we derive 01.7p in the sense of distributions 
on Q. From 02.12p we infer that ip is bounded since W C L°^{VL). As A is continuous, the 
same holds for A( 9 ?), so that X{ip)dtip G L^(0,T;i7). By comparison in 01.71) . Ai? belongs 
to L^(0,T;i7) and 01.71) holds a.e. in Q. It also follows that the normal derivative dn'& 
makes sense in a proper Sobolev space of negative order on the boundary and that the 
integration-by-parts formula holds in a generalized sense. By applying it, one immedi¬ 
ately derives the dynamic boundary condition contained in 01.9p in a generalized sense. 
By comparison in it, dn'd G L^(E) and the boundary condition holds a.e. on S. Finally, 
we remark that 02.10p and the trace theorem imply ■dr £ L°°(0, T; Fr^/^(F)). 

Our hrst result ensures well-posedness with the prescribed regularity, stability and 
continuous dependence on the control variable in suitable topologies. 

Theorem 2.2. Assume 02.3p - 02.51) and 02.91) . Then, for every u G LfifT), problem 
02.14p - 02.17l) has a unique solution (d, dp, C) satisfying 02.10l) - 02.13p . and the estimate 

||'*9||Rl(0,T;R)nL°°(0,T;y) + |||| (0,T;Rr) 

+ ll<dlliyL°°(0,r;R)nRl(0,T;V)nL°°(0,T;W) + ||^||L°°(0,r;R) < Ci (2.18) 

holds true for some constant Ci that depends only on O, T, the structure 02.31) - 02.5p of 
the system, the norms of the initial data associated to 02.9p and an upper bound for ||m|| 2 - 
Moreover, if Ui G i = 1,2, are given and (19*, 6) Ihe corresponding 

solutions, then the estimate 


||'*9l — 'd2||L°o(o,r;R) + 11(1 * '(^l) — (1 * 'd2)||L°°(0,T;y) 

+ ||'*^i,r — 'd2,r||L2(s) + IIt’i — T2\\m(o,T-,H)nL<=^(o,T-,v) 

< C'||ri — R2||l2(s) (2-19) 

holds true where C depends only on O, T and on the structure (1231)-(123]) of the system. 

Remark 2.3. Since W is compactly embedded in (^^(O), the space of continuous functions 
on Q, the regularity 02.12p implies ip G C^{Q) := G°([0, T]; G°(r2)) and the estimate 
ll<d||oo A c||‘d||L°°(o,T;VL)) where c depends only on O. Therefore, we also have 

cp G C^iQ) and \\ip\\co < C 2 (2.20) 


where G 2 is a multiple of the constant Ci of 02.18p . 
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Some further regularity of the solution is stated in the next result of ours. 

Theorem 2.4. Under assumptions fl2.3j) - fl2.5jl and fl2.9jl . the following properties hold 
true. 

i) If in addition Dq G L°°{Q) and u G we also have 

G L^iQ) and ||i9||oo < C 3 (2.21) 

where C 3 is a constant with the same dependencies as Ci and depending on 111901100 o-nd 
on an upper bound for ||m||oo; in addition. 

ii) By also assuming (3°{(po) G we have that f G L°°{Q) and 

llelU^(Q) < (2.22) 

with a constant C 4 that depends on the norm ||/S°(v5o)||oo ns well. 

The well-posedness result for problem fl2.14p - fl2.17p given by Theorem 12.21 allows us 
to introduce the control-to-state mapping S and to address the corresponding control 


problem. We dehne 

X ;= L°°(S) and Vi x ^2 x ^3 where (2.23) 

:= {n G L2(g) ; 1 * n G L“(0, T; 1")}, % := L\T.) (2.24) 

and Vs :=C'°([0,T];i7)nL2(0,T;l/) (2.25) 

S : X —)■ y, M H-)■ S(m) =: ('d, dr, (p) where 

(d, dp, 93,0 is the unique solution to fl2.lUp - fl2.17p corresponding to u. (2.26) 

Next, in order to introduce the control box and the cost functional, we assume that 

Mmin,Mmax ^ L°°(S) Satisfy Wmin < ^max a.e. on S (2.27) 

Ki, e [0,+oo), dg G L^(g) and cpn E (2.28) 

and dehne and 3 according to the Introduction. Namely, we set 

had • ^ ^ • Ilmin ^ H ^ Umax a.e. On Sj’ (2.29) 

3 := ° S : ^ ® where V ^ 1^ is dehned by (2.30) 

’’r. - Y / + Y ^ MT) - (2.31) 


Here is our hrst result on the control problem. 

Theorem 2.5. Assume fl2.3p - fl2.5p and fl2.9p . and let Uad nnd 3 be defined in fl2.23p - 
fl2.3ip . Then, there exists u* G Uad such that 

3{u*) < 3iu) for every u G Uad- (2.32) 


From now on, it is understood that the assumptions fl2.3p - fl2.5p and fl2.9p on the 
structure and on the initial data are satished and that the map S, the cost functionals 
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So and 0 and the control box Uad are defined by fl2.23p - fl2.3ip . Thus, we do not remind 
anything of that in the statements given in the sequel. 

Our next aim is to formulate first order necessary optimality conditions. As Uad is 
convex, the desired necessary condition for optimality is 

(DS(u*),u — u*) > 0 for every u G Uad (2.33) 

provided that the derivative D3{u*) exists in the dual space X* at least in the Gateaux 
sense. Then, the natural approach consists in proving that S is Frechet differentiable at 
u* and applying the chain rule to 3 = So o S. We can properly tackle this project under 
some further assumptions that are satisfied for each of the potentials fll.3p - fll.4l) . We also 
have to require something more on A. Namely, we also suppose that 

D{(3) is an open interval and is a single-valued on D{(3) (2.34) 

/3, TT and A are functions. (2.35) 

In particular, (3° = (3. Furthermore, the inclusion in fl2.15ll reduces to ^ = (3{(p), and 
it is no longer necessary to split the nonlinear contribution to the equation in the form 
^ + 71 {(f). Hence, we set for brevity 


-y := (3 + 71 (2.36) 

and observe that 7 is a function on D{f3). 

As assumptions fl2.3411 - 02.351) force (3{r) to tend to ±cxo as r tends to a finite end¬ 
point of D{(3), if any, we see that combining the further requirement fl2.34p - fl2.35p with 
the boundedness of (p and ^ given by Theorems 12.21 and 12.41 immediately yields 

Corollary 2.6. Under all the assumptions of Theorem fS.fl suppose that 02.341) - 02.35p 
hold, in addition. Then, the component p of the solution (■d, ■dp, 7^,0 satisfies 

P» ^ P ^ P* in Q (2.37) 

for some constants p, ,p* G D{(3) that depend only on G, T, the structure 02.31) - 02.5p and 
O 2 . 34 p - O 2 . 35 p of the system, the norms of the initial data associated to 02.9p . the norms 
ll'dolloo o-nd ||/5(<yCo)lloo) and an upper bound for ||m||oo- 


As we shall see in Section |5l the computation of the Frechet derivative of S leads to 
the linearized problem that we describe at once and that can be stated starting from a 


generic element m G X. Let m G X and h G X be given. We set ('d,'dr,p) '■= S(u). Then 
the linearized problem consists in finding ( 0 , ©p, $) satisfying 

e e H^(0,T;Id)nL^(0,T;V) (2.38) 

0p G 771 ( 0 , T;i7p) and 0p(t) = 0(t)|^ for a.a. 7 G (0, T) (2.39) 

$ G 771 ( 0 , T; 77) n L°°(0, T; V) C 7:^(0, T; W) (2.40) 
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and solving the following problem 


dtQv+ / V0-Vp+ / A((/9)c?f<h T + / 


+ T 


dtQr vy + 01 


©r Pr — Oi 


mhvr 


Jr Jr Jr 

for every (p,Pr) ^ V x Vp such that Pr = P|p and a.e. in (0,T) 
(9t<h — crA$ + 7 '(^) $ = 'dA'(<^) $ + A(<^) 0 a.e. in Q 
dn^ = 0 a.e. on S 
0(0) = 0 and <h(0) = 0 a.e. in hi. 


(2.41) 

(2.42) 

(2.43) 

(2.44) 


Proposition 2.7. Letu G X and {'d,'dr,lf) = S(p). Then, for every h G X, there exists a 
unique triplet (0, ©p, $) satisfying fl2.38p - fl2.40p and solving the linearized problem fl2.4ip - 
fl2.44p . Moreover, the inequality 


||(0,0p,$)||y <C5||h||x (2.45) 

holds true with a constant that depend only on hi, T, the structure fl2.3p - fl2.5p and 
fl2.34p - fl2.35l) of the system, the norms of the initial data associated to fl2.9l) . and the 
norms ||p|1ooj ||'*^o||oo and ||/S(v5o)||oo- la particular, the linear map T> ■. h ^ (0, ©r, $) is 
continuous from X to V. 


Namely, we shall prove in Section [5] that the map T> G X(X, V) introduced in the last 
statement exactly provides the Frechet derivative DS(u) of S at u. Once this is done, 
we may use the chain rule with u := u* to prove that the necessary condition fl2.33p for 
optimality takes the form 

Ki [ if)* -’§0)0) + K 2 [ (v 3*(T) - v9n)*h(T) > 0 for any M G Had (2.46) 

JQ Jn 

where (i?*, 'df, ip*) = S(p*) and, for any given u G Uad, the triplet (0, ©r, $) is the solution 
to the linearized problem corresponding to h = u — u*. 

The hnal step then consists in eliminating the pair (0,$) from fl2.46p . This will be 
done by introducing a triplet {p,pr,(l) that fulhls the regularity requirements 

p,q(E H\0,T-,H)nL°°{0,T-,V), q e L\0,T-,W) (2.47) 

Pr & Hp) and pp{t) = p{t)\^ for a.a. f G (0, T) (2.48) 

and solves the following adjoint system: 


- / dtpv+ / Vp-Vp- / X{Lp*)qv - T / dtppvp + a / prPr 
Jn Jn Jn Jr Jr 


= Ki {- d * - i9q)p 
Jn 

for every v Vp := P|p, and a.e. in (0,T) 

- [ dtqv + a j Vg-Vp+ f -'d*X'{ip*))qv = j X{^p*)dtP'^ 

Jn Jn Jn Jn 

for every n G P and a.e. in (0,T) 
p{T) = 0 and q{T) = K 2 {(p*{T) — (p^) a-e. in 12. 


(2.49) 


(2.50) 

(2.51) 
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Clearly, fl2.49jl - fl2.50j) are the variational formulation of a boundary value problem. 
Namely, p and q solve two backward parabolic equations complemented by a dynamic 
boundary condition for p and the homogeneous Neumann boundary condition for q. How¬ 
ever, it is more convenient to keep the problem in that form. 

Theorem 2.8. Let u* and (i?*, "dp, 9?*) = S{u*) be an optimal control and the correspond¬ 
ing state. Then the adjoint problem fl2.49p - fl2.51l) has a unique solution {p,pr, q) satisfying 
the regularity conditions fl2.47l) - fl2.48p . 

We recall that, if K is a closed interval and yo G K, the normal cone to K at yo is the 
set of z G M such that z{y — 1/0) < 0 for every y E K. Here is our last result. 

Theorem 2.9. Let u* be an optimal control and ('d*,'dp,(p*) = S(u*) denote the associate 
state. Moreover, let {p,pr,q) be the unique solution to the adjoint problem fl2.49l) - fl2.5ip 
given by Theorem A2.8\. Then, for a.a. {t,x) G S, we have 

m{x)pr(t, x) {y - u*(t, x)) > 0 for every y G [Mmin(t, 2:), Mmax(C a:)] (2.52) 

that is, —m{x)prit,x) belongs to the normal cone to [umwit, x), Rma.x(h x)] atu*{t,x). In 
particular, we have 


U Rmax) R Rmin and U G [itmin) Rmax] 
a.e. in the subsets ofE where mpr < 0, mpr > 0, mpr = 0, respectively. 

In performing our a priori estimates in the remainder of the paper, we often use the 
Holder inequality (with the standard notation p' for the conjugate exponent of p), its 
consequences and the elementary Young inequalities 

ab < uj -|- (1 — w) gP) < ^^2 _|- pi 

46 

for every a,b >0, ca G (0,1) and 5 > 0 (2.53) 

as well as the continuous (in fact compact) embedding V C L^(H). Moreover, in order 
to avoid a boring notation, we follow a general rule to denote constants. The small-case 
symbol c stands for different constants which depend only on H, on the hnal time T, the 
shape of the nonlinearities and on the constants and the norms of the functions involved in 
the assumptions of our statements. A small-case symbol with a subscript like cs indicates 
that the constant might depend on the parameter 6 , in addition. Hence, the meaning 
of c and cs might change from line to line and even in the same chain of equalities or 
inequalities. On the contrary, we mark precise constants which we can refer to by using 
different symbols, e.g., capital letters. 


3 The state system 

This section is devoted to the proofs of Theorems 12.21 and 12.41 As far as the former 
is concerned, we start proving its second part, i.e., the continuous dependence formula 
fl2.19p . From this we derive uniqueness as well. 
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Continuous dependence and uniqueness. We first derive an identity that is satisfied 
by any solution. By integrating fl2.14p with respect to time, we obtain 


'dv+ V(l*'d)-Vn + r / ■drur + « / * {'dr —'rnu))vr 

Jn Jq Jr Jr 

= [ 'dov + T [ (i9o|r)Pr- [ P(<^) - ^(</?o))p 
Jn Jr Jn 

for every (n,nr) & V x Vr such that ur = P|p and a.e. in (0,T). (3.1) 


Now, we pick two elements Ui G i = 1, 2, and consider two corresponding solutions 

{'di,'di^r,^i,^i)- We write fl3.ll) for both controls and solutions and test the difference by 
choosing v = -d := 'di—'d 2 and vr = 'dr ■= Then, we integrate over (0, t), where 

t G (0, T) is arbitrary. At the same time, we write fl2.15l) for both solutions, multiply the 
difference by (f := (fi — (f 2 and integrate over Qt- Finally, we add the equalities we have 
obtained to each other. By also setting u := Ui — U 2 and ^ ■= — ^2 for brevity, we 

infer that 


/ + l /|V(l*tf)(t)P + fo + | /l(l*tfr)(t)P 

j Qt St 

+ R / + ^ [ |Vv?P+ f 

^ Jn JQt JQt 

= a f m{l*u)'dr- f (X{(pi) - y^{cp 2 ))'d 
J^t Jqi 

- f {ti{lpi) - ti{lp2))lp + ( {'di\{ipi) -'d2X{(p2))^- (3.2) 

^ Qt J Qt 


All the terms on the left-hand side are nonnegative, including the last one since jS is 
monotone. We estimate each term on the right-hand side, separately. In the sequel, 6 is 
a positive parameter. We have 


a / m{l * u)'dr =—<y / mu{l*'dr) + <y / m{l * u){t) {1 *'dr){t) 


'St 


'St 


<c / |l*'drP + c / + 5 / 1(1 * ■dr)(^)P + C5 / \{^*u){t)\‘^ (3.3) 


'St 


'St 


and the last integral is bounded by C5 ||m||p 2 (s) to the Holder inequality. Next, owing 
to the boundedness of and Lp 2 ensured by Remark 12.31 and to the regularity of A on 
bounded intervals, we infer that 


'Qt 


(X{^Pl) - Xf2))^ < C / |(p| |l9| < (5 


'Qt 


'Qt 





2 


and the third integral on the right-hand side of fl3.2p can be treated in a similar way. 
Finally, we deal with the last term. As A is Lipschitz continuous (see fl2.5lB and (p 2 is 
bounded, we have 


'Qt 


(i9iA(</?i) - A(<^2))</?= / 'di{X{(pi) - X{(P2))^ + 'd\{f2)'P 


' Qt 


' Qt 


<c |i9i| + c / |i9| |<^| < c / |i9i||<^p + (5 


+ Cs 


'Qt 


'Qt 


'Qt 


^ Qt 


^ Qt 
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On the other hand, by combining the Holder ineqnality, the Sobolev ineqnality ||n||4 < 
c||n||y for every v E V, and the regnlarity -di E L°°(0,T; V), we obtain 

[ [)\Ms)\UMs)\UMs)hds 

JQt Jo 

< C [ \Ms)\\v\Ms)\\ 2 ds < 6 [ |V(^P + C5 /" 

^0 ^ Qt ^ Qt 

At this point, we collect all the estimates we have derived, choose 6 small enough and 
apply the Gronwall lemma. Thus, we obtain fl2.19p and uniqueness easily follows. Indeed, 
taking ui = U 2 in 02.191) immediately yields -di = '^ 2 , 'dip = 'd 2 ,r and (pi = ip 2 . By 
comparison in 02.151) . we also have = ^ 2 - Gl 


In order to complete the proof of Theorem 12.21 we have to show the existence of a 
solution and to establish estimate 02.181) . To this end, we hrst replace /9, A and A by the 
smooth approximation of them fde, Ag and A^ we introduce below, where e is a positive 
parameter, say e E (0,1). By doing that, we obtain the approximating problem of hnding 
a quadruplet {'de,'de,r,^e,ie) satisfying regularity requirements of type 02.1Up - 02.13p and 
solving 

I dt'deV+ I Vde-Vn+ f dtX{^e)v + T I dt-de^Wr + a f {'d^^r - mu) vr = 0 


Jn Jn Jo. Jv Jv 

for every {v,vr) E V x Vr such that ur = and a.e. in (0,T) (3.4) 

dt^e - + Ce + = ’deXXe) and = XXe) a.e. in Q (3.5) 

dn<^e = 0 a.e. on E (3.6) 

'de(O) = 'do and <yC£(0) = (po a.e. in H. (3.7) 

The regularity we require for the solution is the following 

E H\0,T;H)nL^{0,T;V) (3.8) 

■dep G 77^(0, T; ihp) and 'dep(f) ='de(f)|p for a.a. f G (0, T) (3.9) 

if, E H\0, T; H) n L'^(0, T; V) C L2(0, T; W). (3.10) 


The lower level fl3.10p with respect to (12.12p has been chosen for convenience. However, 
the solution we hnd also satishes fl2.12p . as it will be clear from the proof. In the above 
equations, X is the Yosida regularization of (d at level £ (see, e.g., [3 p. 28]). It is well 
known that X is maximal monotone, single-valued and Lipschitz continuous. We also 
introduce the function (de dehned by 

pr 

jde{r)'.= / jde{s)ds for r G M (3.11) 

Jo 


and recall that 


IXiXl A |/d°(^")| and /Se(r) /S°(r) for r E D{(d) 
0 < (deiX A (d{r) for every r G M. 


(3.12) 

(3.13) 
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For fl3.12j) see, e.g., O Prop. 2.6, p. 28], while fl3.13p follows from fl3.12p and /S£(0) = 0 
(cf. fl2.4p ). Furthermore, as well as its derivative Ag are dehned by 

A£(r) := X{r)({sr) and Ae(r) := — Xi;{r) for r G M, where 

CL! 

C e satishes C(^) = 1 for |r| < 1 and ({r) = 0 for |r| > 2. (3-14) 

Notice that both A^ and A^ are bounded and Lipschitz continuous, and we set 

Ae:=sup|A£| and Ae:=sup|A£|. (3.15) 

Our project is the following: i) we prove that problem fl3.4p - fl3.7p has at least a 
solution by a hxed point argument; ii) using compactness and monotonicity methods we 
show that its solution tends to a solution of problem fl2.14^ - 02.17p as £ 0, at least for a 

subsequence. We need two lemmas. 

Lemma 3.1. Let ■dg G L^(0,T;i7). Then, there exists a unique satisfying fl3.10p . 
fl3.5p - fl3.6l) and the second Cauchy condition in fl3.7p . Moreover, the estimate 

ll9^e||Hl(0,T;H)nL°°(0,T;U)nL2(0,T;IT) < ^^(l + \\'de\\L^{0,T-,H)) (3.16) 

holds true with a constant that depends on the structure fl2.3p - fl2.5lh the norms involved 
in (O and e, but it is independent of 'de- 


Proof. We are dealing with a standard semilinear parabolic problem that has a unique 
solution with the required regularity. We just derive estimate fl3.16p and control the 
dependence of constants. We multiply equation fl3.5l) by dt<Pe and add the same integral 
to both sides, for convenience. We have 


a 


i^.r + 77 \^Ts{t)\^+ ^siPeit)) 


' Qt 


a 

2 


J Qi ^ J ^ J ^ 

|V<PoP+ / Pe{Po)- I '^(Te)dt(Pe+ [ KiPe) dt(Pe + [ We 

jQt 


(3.17) 


' Qt 


' Qt 


where is given by fl3.1ip . We recall fl3.13p and fl2.4l) . which imply that jdeWo) < 
PWo) < f3°Wo)To a.e. in hi. Thus fl2.9p yields ||;d£(<po)||i < c- Furthermore, notice that 
|A£((pe)| < Aj (see fl3.15p j. On the other hand, vr has a linear growth, so that 


7r{ipe)dtiPe < 


IdtPe]"^ d-C / + C 


^ Qt 


^ Qt 


^ Qt 


and the last integral, which coincides with the last term on the right-hand side of fl3.17l) . 
can be treated as follows: for every s G [0,f], we have 

(p£(s)=v9o+ / dt^PeW) ds', whence 
Jo 

WeW)\^ <2Wo?+ ‘^ [ dt(Pe{s')ds' <c + c f \dtipWs')\'^ dsC 

Jo Jo 

It turns out that 

Wsit)]"^ < C + c [ IdtPel'^ and [ We\‘^<c + c[ ([ ds . (3-18) 

./Qj In. .In ^.In ' 


'O 


' Qt 


'0 ^JQs 
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The second inequality in fl3.18jl implies that its left-hand side, i.e., the integral we are 
dealing with, can be handled by the Gronwall lemma and we conclude that 

\\Ve\\H^{0,T-,H) + II V<y5£||2,oo(o,T;R) < C + c(l -|- ||'d£||L2(o,T;R) 

where c depends only on the structure and the initial datum ipQ. From this estimate and 
the hrst inequality in 03.181) it follows that 

ll9^£||Rl(0,T;R)nL°°(0,T;T/) < C + c(l -|- Ag) ||||L 2(0,r;R) 

with a similar new constant c. Now, for a.a. t G (0,T), we write 03.5p at time t, multiply 
by ^^{t) and integrate over fl. We obtain a.e. in (0,T) 

f = / {-dt^e - K{^e)) is 

Ju Ju Jo. 

— of + C + ^ + + \'&eK{^s)\‘^) 

^ Jq Jn 

whence immediately for a.a. t G (0,T) 

IIC^WIIr < c(l + \\dt(pe{t)\\H + \Mt)\\H + |li9£(t)A£((p£(f))||/f). (3.19) 

In view of the regularity of (pe already achieved and the one of recalling 03.151) we 
deduce that 

ll^£llL2(o,r;iy) < c -1- c (1 -|- Ae) ||'d£||2,2(o,r;iy) 

where c has the same dependencies as above. As the estimate of the norm of in 
L^(0,T;hF) follows by comparison in 03.5p and elliptic regularity, the proof of 03.16p is 
complete. □ 

Lemma 3.2. Let pe ^ H). Then, there exists a unique pair ('de,'d£,r) satisfying 

03.8p - 03.91) and the first initial condition in 03.71) and solving the variational equation 03.4p . 
Moreover, the estimates 


\mLHQ)<Re (3.20) 

||'*^£||Rl(0,T;R)nL°°(0,T;V) + ||||(0,T;Rr) ^ -^e(l + IIII(0,T;R)) (3.21) 

hold true with constant and that depend on the structure (CT-ira. the norms 
involved in 02.91) . the norm ||m ||2 and e, but they are independent of pe. 


Proof. We set 

V := {(n,nr) E V x Hr : vr = n|p}, H := H x Hr 
and endow these spaces with the scalar products dehned on and by 

((ta, tcr), (L’,L’r))v := (fi7w-Vv + wv)+ wrVr 

Jn Jr 

{{w,wr),iv,vr))j^:= / wv + r wrVr 


(3.22) 


(3.23) 

(3.24) 
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respectively. Then, we obtain two Hilbert spaces and V is continuously and densely 
embedded in !K, so that we can construct the Hilbert triplet (V, Tt, V*) in the usual way. 
Moreover, we dehne the continuous bilinear form a on by the formula 

a((ta, tcr), (p, Pp)) := / Vta ■ Vn + a / tcrur (3.25) 

Jn Jr 

and consider the operator A G /C('V, V*) associated to a. As a clearly satishes 

a((n,nr), (u,nr)) + (l + 7)||(wAr)||M > ll(wAr)llv for every {v,vr) G V 

the general theory (see, e.g., m ) ensures that, for every F G L^(0,T;'V*) and Uq G IK, 
there exists a unique U satisfying 

U G H\0,T;V*)nL‘^{0,T;V) C ^“([O, T]; A) 

U'{t) + A U{t) = F{t) a.e. in (0, T) and f/(0) = Uq 

and that U G H^{0,T;A) D L°°{0,T;V) whenever F G L^(0,T;1K) and Uq G V. By 
accounting for our assumption on ip^ and fl2.9p . we choose 

F = {-dtX{pe),Oimu) G L^(0,T;A:) and Uq = {'iJo,'do\X ^ 

and obtain a unique pair ('dg, 'de^r) satisfying fl3.8l) - fl3.9p and the hrst initial condition fl3.7p 
and solving the variational equation fl3.4p . Let us now prove estimates fl3.20p - fl3.2ip . We 
observe that the analog of fl3.1l) obtained by replacing A by holds for (-dg, So, 

we take ('d£,'de,r) as test pair and have 

f + l /|V(l»il.)(t)r + r / /|(l»tf..r)wr 

JQt ^ Jn JT,t ^ dr 

= a[ m{l*u)iJe,r+ [ A'^e + rf (i9o|r)4,r + / (^(t’o) - • (3.26) 

dst J Qt J "Tt J Qt 

All the integrals on the left-hand side are nonnegative and the hrst term on the right-hand 
side can be estimated as we did in fl3.3p . namely 


af m{l *u)'ds,r < c f \l * + c f \u\‘^ + 6 f \{1 * 'ds,r)it)\‘^ + cs [ * u)(t)\‘^ 

J^t J'Tt Jr,t Jr Jr 

and the last integral is bounded by C5 ||m||^ 2 (s) due to the Holder inequality. The remaining 
terms on the right-hand side of fl3.26l) are treated in the usual way with the help of the 
Holder and Young inequalities; just for the last term we point out that 


(X(<^o) - XXe))'Je < 2As / <S + C^A^ 


^ Qt 


^ Qt 


' Qt 


thanks to fl3.15l) . Thus, by choosing 5 small enough and applying the Gronwall lemma, 
it is straightforward to obtain the desired estimate fl3.20l) . Let us now prove fl3.21l) . The 
rigorous argument could rely on testing fl3.4p by a V-valued approximation of {dt'dei dt'deX- 
However, we prefer to avoid such a detail and formally test the equation by (St'de, dt'de,r)^ 
directly. We have 


a 


' Qt 


WeV + r: / |Vi9£(t)|^ + r / + ^ / |4,rWI 


'St 


= - / Xe{ipe)dtPedt'Je + O: / mU + X / |Vl9oP + 


a 


' Qt 


'St 


l^oh 























Colli 


Gilardi 


Marinoschi 


15 


As |Ae(</?£)| < Ae (see (13.15^ 1 and "do satisfies fl2.9p . we immediately derive fl3.2ip owing 
to the Holder and Young inequalities. □ 

At this point, we are ready to complete the proof of Theorem 12.21 by following our 
project sketched above. 

Existence of the approximating solution. As said before, we are going to use a 
hxed point argument. We often avoid stressing the dependence on £, which is hxed. We 
consider the closed ball of L‘^{Q) 

‘B-={ve L'^iQ) : |ln||2 < Re} (3.27) 

where is given by Lemma |3Y] (see (13.201) 1 and dehne the map T : (B —)■ L^{Q) by the 
following steps: i) for G !B we apply Lemma 13.11 where is replaced by d, hnd the 
solution and term it £(1?); ii) by starting from such a we apply Lemma [3.21 hnd 
the solution ('de,'de,r) and set T('d) := dg. By construction, it turns out that T('d) G ‘B: 
indeed, the constant Re in (13.201) is independent of \\pe\\m{o,T-,H)- Moreover, with the 
above notation, we deduce from Lemmas 13.11 and (cf. (I3.16P and (13.211) 1 

11 1 1 (0,T;R)nL°° (0,T;y)nL2 {0,T;W) < Ce{l + Re) (3.28) 

|l'*^£:||Rl(0,T;R)nL°°(0,r;V) + ||'*^£:,r ||(0,T;Rr) ^ + Ce{^ + -^e)) • (3.29) 

Therefore, T(!B) is relatively compact by the Aubin-Lions lemma (see, e.g., [25l Thm. 5.1, 
p. 58]). Now, we verify that T is continuous. So, we assume that G (B and that 
converges to d in L^(Q), and we prove that := T('d„) converges to T('d) in L‘^{Q). 
We set for convenience := £('d„) and := d^ip. As G (B, estimates (13.281) - 
(I3.29P hold for and d„_r- Hence, for a subsequence, d„ and converge to 

some p, d and dr in the correponding weak or weak star topologies. Clearly, dr = d|p. 
Moreover, pn and d„ converge to p and d strongly in L‘^{Q) by the Aubin-Lions lemma. 
This implies that f{pn) converges to f{p) strongly in LF‘{Q) for any Lipschitz continuous 
function / : M —)■ M, and this is the case if / is either Ag, or Ag, or de, or tt. As a 
consequence, dtXe{pn) = Xe{pn)dtPn and d„Ag(<yC„) converge to Xe{p)dtp and dAg(:p) at 
least weakly in L^{Q). Therefore, the quadruplet (d,dr,V5,0) where ^ := Pe{^), solves 
the integrated version of the variational formulation of the approximating problem (13.41) - 
(13.7D with smooth test functions, with dg replaced by d in (13.51) . This easily implies that 
d = T(d). As the same argument holds for any subsequence extracted from {d„}, the 
continuity we have claimed is proved. Therefore, we can apply the Schauder hxed point 
theorem and conclude that T has a hxed point dg. As dg belongs to T((B), it satishes (13.81) . 
Moreover, dg,r := dg|p belongs to Hr), the function pe := £(dg) satishes (I3.10p 

and the quadruplet (dg, dg,r, </?£, ^e), where ^g := (3e{pe), solves (I3.4p - (l3.7p . □ 

The last step consists in letting e tend to zero and getting a solution to (I2.14p - (l2.17p . 
To this aim, we derive a priori estimates on the approximating solution that are uniform 
with respect to e G (0,1). As such estimates are conserved in the limit as e n 0, inequality 
(I2.19P is established as a consequence. 

First a priori estimate. We choose n = dg and ur = '&£,r in (13.41) and integrate 
over (0,f). At the same time, we multiply (13.5p by dtPe and integrate over Qt. Then, we 
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add the equalities obtained in this way and observe that the terms involving the nonlinear 
function cancel out. Hence, we have 

^ [ l^£(^)r+ [ + 7^ f + a [ 

^ Jn JQt ^ Ast 

+ [ \dt^e\^ + ^ [ |V^.(t)r+ [ dei^eit)) 

= 2 / l'*^oP+2 f l'*^o|rP + ^ [ |Vv9oP+ f 

+ a mui!}s,r - / Ti{^e)dt^e ■ 

JT,t JQt 

As fl3.13p holds, tt is Lipschitz continuous, m and u are bounded, and the data satisfy fl2.9l) . 
we easily deduce that 

||'*^e||L°°(0,T;H)nL2(o,r;U) + H'^eT|U°°(0,T;Hr) + ||||Hi(0,T;iy)nL°°(0,T;U) < C. (3.30) 

Since |A£(r)| < |A(r)| + c(l + |r|) < c(l + |r|) (by fl3.14p and fl2.5p l and V C h®(f2), it 
follows that 


\\K{.¥’e)\\L^{o,T-,Lf^(n)) ^ c and \\'d£X£{(p£)\\L^{o,T-,L3{n)) < c. (3.31) 

Therefore, by applying fl3.19p . we deduce that 

ll'Ce||L2(0,T;H) < C (l + ||V?e||(0,T;H) + \\‘^eK{^e)\\{0,T;H)) < C. (3.32) 

By comparison in 03.51) . we derive an estimate of Aip^ in L‘^{0,T] H), whence 

|Iv^£||l2(o,T;IT) < c (3.33) 

by 03.op and elliptic regularity. 

Second a priori estimate. The estimate we need next is the following 

||^£||Hl(0,T;H)nL°°(0,T;y) + || I?£,r || Rl (0,r;Hr) 

+ ll<^£|l ITL°“(0,T;H)nHl(0,T;y) < C . (3.34) 

A rigorous proof is given in the Appendix. Here, we proceed formally. We take v = dt'&£ 
and vr = as test functions in 03.4p and integrate over (0,t). At the same time, we 

multiply the equation obtained by differentiating 03.5p with respect to time by dt^Pe and 
integrate over Qt. Then, we add the equalities just derived to each other. Since the terms 
involving the product dj'de dtip£ cancel out, we have 

/ mf + U\^MtP + r [ \m,rr-+'^ 

+ l [ \\/dtp>£\‘^+ [ l3'^{^Q\dtP>£\‘^ 

^ Jn JQt JQt 

= \ j iV-doP + l f l^o|rP + ^ [ \dt^e{0)\‘^+ j ^£K{^£)\dtip,\^ 

+ a mu dt'&£^T — 

JT.t 


Ti'icps) 


' Qt 
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All the terms on the left-hand side are nonnegative, in particular for the last one we use 
the monotonicity of (3^. Moreover, the hrst two terms on the right-hand side are hnite by 
02.91) and the last two integrals can be treated in an obvious way by also taking the Lip- 
schitz continuity of tt into account. In order to control the term involving A', we observe 
that 03.14P easily implies that |A^(r)| < c for all r G M and e G (0,1). We also recall that 
V C L'^{Q) and term C the norm of the embedding. Therefore, by owing to the Holder 
and Young inequalities, we obtain 

[ 'deKiV>e)\dt^e\'^ < C [ ||l9£(s) jU || 5 jV 3 £(s) II 4 ||ai99<,(s) jja ds 
JQt Jo 

Jo + \\Ms)\\l\\dt^e{s)\\lds 

/ \\^^t^e{s)\\Hds + ^ \\dtiPe{s)\\H ds + C \\lJe{s)\\l \\dtiPeis)\\l ds 

and we remark at once that s i-A- H'd^ is bounded in L^(0,T) by fl3.30p . Finally, we 
estimate dt(pe{0) in H. We formally have from fl3.5l) 

dt<Pe{0) = aA^po - (JsM - 7r((^o) + ^oA£((^o) a.e. in n 

and we can owe to fl2.9p . fl3.12l) . the inequality |A£(r)| < |A(r)| -|- c(l -|- |r|) < c(l -|- |r|) 
for every r G M, and account for the continuous embedding V C L'^{Q) once more. Then, 
fl3.34p follows by the Gronwall lemma. 

Third a priori estimate. The Holder inequality, the continuous embedding V C L^{Q) 
and fl3.34l) imply 

||'deA£((^e)||Loo(o,T;R) < C ||•de||Loo(0,T;L4(^2)) (l + ||‘^e||L°°(0,T;L4(f2))) < C. (3.35) 

On the other hand, fl3.19p holds for the approximating solution. Thus, is bounded in 
L°°{0,T-, H) and a bound for in L^{0,T-, H) follows by comparison in fl3.5l) . Hence, 
by also using fl3.6l) and the elliptic regularity theory, we have 

\\^£\\l°°{0,T-,W) + ||'Ce||L°“(0,T;R) < C- (3.36) 


Conclusion of the proof. By standard weak and weak star compactness results, we 


have for a subsequence 



weakly star in H^{0, T; H) D L°°(0, T; V) 

(3.37) 


weakly in 77^(0, T;i7r) 

(3.38) 


weakly star in IFi’°°(0, T; 77) n77i(0,T;I/) nL°°(0,T; W) 

(3.39) 


weakly in L°°(0, T; 77) 

(3.40) 


and iJrit) = 'd(t)|p for a.a. t G (0,T). By owing to the compact emedding W C C°(f2) 
and to m Sect. 8, Cor. 4], we can also assume that the selected subsequence satishes 


ipe ^ ^ uniformly in Q 


( 3 . 41 ) 
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SO that converges to A(<yc) uniformly in Q since Ae(r) converges to A(r) uniformly 

on every bounded interval. We deduce that X^{ips)dt(p£ and converge to X{ip)dt(p 

and to "dX^ip) at least weakly in L‘^{Q). Finally, we have ^ G P{(p) a.e. in Q by applying, 
e-g-, [D Prop. 2.2, p. 38]. Therefore, we can pass to the limit in the integrated version of 
problem fl3.411 - 03.71) written with time dependent test functions and easily conclude that 
(•d, dr, V?, 0 is a solution to problem 02.14^ - 02.171) . □ 

Now, we prove Theorem 12.41 For the claim i), we hrst consider the Cauchy problem 
for the linear variational equation 


dt'&v + / Vd ■ Vn + T dt'&r vr + a / dr ur = / i/jv + / i/JrVr 


for every {v,vr) G C x W such that up = P|p and a.e. in (0,T) 


(3.42) 


where and ipr are given, and prove the following 

Proposition 3.3. Assume that ip G L°°{0,T; H), ipr G L°°(S) and do G P nL°°(r2) and 
that the corresponding norms are estimated by some constant C. Moreover, assume that 
(d,dr) enjoys the regularity 02.101) - 02.lip and satisfies the variational equation 03.421) and 
the initial condition d(0) = do- Then, d is bounded and the estimate 

||d||oo<C 

holds true with a constant C that depends only on hi, T, r, a and C. 


Proof. This is a regularity result for a linear problem. Hence, it can be established by 
considering the following cases: 

a) Ip = 0 and do = 0; b) ipr = 0 . 

Let us consider the hrst one. For any integer n > 0 and p G (2, +cxo) we dehne 

r ~ /■I'’! -1 

T„(r) := ntanh — and Tnp{r) := J (T„(s))^ ds for r G M. 

Now, we set d^ := T„(d) and d„,r := Fn(dr), and we test 03.421) by n = |d„|^“Wignd 
and Ur = |dn,r|^~^signdp, where the sign function is extended by sign(O) = 0. Then, we 
integrate over (0,T). As signd^^r = sign dp, we obtain 


T„p(d(T)) + r / f„p(dp(T)) + (p-l) / |d„r2T'(d„)|Vd|2 + a / |dp||d„,pri 


'Q 


= / V^rldn.pl^ ^ sign dr 


All the terms on the left-hand side are nonnegative. Hence, by observing that |r| > |Tn(p)| 
for every r, and applying the Young inequality, we have 


« / |d„,pr<a / |dp||d„,pri< / \iPr\K,rr^ 






















Colli 


Gilardi 


Marinoschi 


19 


whence immediately 


a 

P 


|i9n,rr < / \i^l 

p 


or 


a 


i/p 


II# 


n,T\\p 


< a 


-w 




rllp' 


By noting that |T„(r)| /> |r| for every r G M as n ^ oo, we can let hrst n and then p tend 
to infinity and dednce that #r is bonnded. Namely 


||#r|loo < a llV'rlU < C/a. 


(3.43) 


Hence, we can apply 1221 Thm. 7.1] with q = 2 and r = oo and conclude. 

In case b) we adapt the proof of [22l Thm. 7.1] (still with q = 2 and r = oo) to our 
situation. For k > max{l, C} we set := (# — k)~^ and := ('dr — k)^ and take v = 
in fl3.42p . By simply writing 'dr = ('dr — k) + k and observing that 'd^(O) = 0 since 'do < C, 
we obtain 




whence also 


|4(t)P+ / |V# 

jQt 


k\2 


T CK 


|4|2 + A;a 


dj; = 




'St 


'St 


'Qt 


1 

2 


1^ + / |vr|"< / 


'Qt 


' Qt 


This corresponds to formula [2^ (7-6), p. 183] and the whole argument can be performed 
in the same way till formula [221 (7.14), p. 186]. Then, just one modihcation is needed. 
Namely, we account for m Rem. 6.2, p. 103] in applying [221 Thm. 6.1, p. 102] since no 
upper bound for 'dr is known now. This leads to the desired estimate from above d <C. 
The corresponding estimate from below is obtained by applying the former to —'d. □ 


Now, we apply the above result by observing that the pair ('d, 'dr) we are interested in 
satishes fl3.42p with := X{p)dtq^ and '^r := iRM, and notice that these functions belong 
to L°°{0,T] H) and to L°°(S), respectively. Moreover, the corresponding norm of ijj has 
been already estimated by fl2.18p . while an upper bound of the norm of 'ipr is supposed 
to be given in the statement. This yields the claim i) of Theorem 12.41 

The next step should be the proof of ii) of Theorem 12.41 To this aim, we just refer 
to m Thm. 2.2, Hi)]. In fact, the proof given there shows that the component of a 
pair [ip, satisfying fl2.12^ - 02.13p and solving the homogeneous Neumann problem for 
the equations 


dt<p — aAip + ^ + vr((p) = / and ^ G /S((p) a.e. in Q 

is bounded whenever / G L°°{Q). Such a statement is proved just with / = 'd in the 
quoted paper by knowing that'd is bounded. However, the same proof is valid with any 
bounded /. Here, we have / = 'dA((p), and both'd and A((p) are bounded. 


4 Existence of an optimal control 


We prove Theorem 12.51 by the direct method. Since Uad is nonempty, we can take a 
minimizing sequence {un} for the optimization problem and, for any n, we can consider the 
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corresponding solution to problem fl2.14p - fl2.17p . Then, {un} is bounded 

in and estimate fl2.18p holds for Therefore, we have 


Un u* 

^n 

■dnT —t 'dp 

‘^n 


weakly star in L°°{Q) 

weakly star in H) nL°°(0,T;l/) 

weakly star in i7^(0, T; Hr) 

weakly star in 1T^’°°(0, T; H) n H\0, T; V) n L°°(0, T; W) 
weakly star in L°°(0, T ; H) 


at least for a subsequence. Then, u* G Uad since Uad is closed in X. Moreover, the 
initial conditions for 'd* and (p* are satished and 'dp = "d^jp a.e. in (0,T). Thus, we can 
easily conclude by standard argument. Indeed, converges to Lf* uniformly in Q 

due to the compact embedding W C (see, e.g., 1211 Sect. 8, Cor. 4]), whence 

n{ipn) and converge to and A(<y5*) in the same topology. We also deduce 

that X{'dn)dtp>n and 'dnX{p>n) converge to \{'d*)dtp>* and 'd*A(<p*) at least weakly in L‘^{Q), 
and that G (note that —)■ strongly in L‘^{Q) and see, e.g., [H Prop. 2.1, 

p. 29]). Hence, ('d*,'dp, (p*, ^*) satishes the variational formulation in the integral form of 
problem fl2.14l) - fl2.17l) corresponding to u*. Therefore 

3{u*)=3o{^*,'&l,(p*)<liniinf3o{'&n,'&n,r,^n)=^i^3{un)= inf 3{u). 

n—^oo n—^oo uGUad 


5 The control-to-state mapping 

As sketched in Sections O the main point is the Frechet differentiability of the control- 
to-state mapping S. This involves the linearized problem fl2.4ip - fl2.44p and we hrst prove 
Proposition 12.71 i.e., well-posedness for the linearized problem and the continuous depen¬ 
dence of its solution on the parameter h. It is understood that all the assumptions of 
Theorem 12.41 as well as fl2.34p - fl2.36p are in force. 

Well-posedness. We aim to apply a contraction argument. To this end, we observe 
that all the coefficients that enter the equations but c := X'(jp)dt'^ are bounded thanks to 
Corollary 12.61 and that the possibly unbounded coefhcient c belongs to L°°{0,T-, H). We 
dehne the maps Ti, T2 and T in a proper functional framework as follows. For 0 G L‘^{Q), 
we consider the problem for $ given by fl2.42l) - fl2.43l) and the second condition in fl2.44p . 
where 0 in fl2.42p is replaced by 0. We obtain a linear parabolic problem which has a 
unique solution $ satisfying fl2.40p . We set Ti(0) := $. By doing that, we obtain a 
map Ti : L‘^{Q) —)■ H). Now, we £x $ G H) and consider the problem 

for (0, 0r) given by the variational equation fl2.41l) and the initial condition 0(0) = 0. 
Such a problem has a unique solution (0, 0r) satisfying fl2.38l) - fl2.39l) . as one can see by 
arguing as in the proof of Lemma [3.21 and we set T2(‘h) := 0. In this way, we obtain a 
map T2 : H^{0, T; H) ^ L'^iQ)- We set T := T2 o Ti and prove that some iterated T™' of 
T is a contraction in L‘^{Q) by deriving some estimates involving Ti and T2, separately. 
For 0 G L^{Q) we write fl2.42p with 0 in place of 0 and multiply it by <h. Then, we 
integrate over Qt and obtain 

f [ ia(»5)i lei i«>|, 

J ^ J Qt J Qt J Qt 
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Hence, with the help of the Young inequality and the Gronwall lemma, we deduce the a 
priori estimate 


c ||0||l2(qo for every t e [0,T] 

where the constant C we have marked with a capital letter for a future reference does 
depend neither on t nor on 0. As 5”i is linear, this means that 

||3^l(0l) — 3^l(02)||L°°(O,t;R)nL2(o,t;V) < c ||01 — 02 ||l2(Qj) for f G [O, T] (5.1) 

for every 0i, 02 G L'^{Q). Now, for given G H), i = 1,2, we consider the 

problems corresponding to the dehnition of 0j = 3^2 (®i) and take the difference. By 
setting $ := $1 — $2 and analogously dehning 0 and 0r for brevity, we see that fl2.4ip 
holds with h replaced by 0. We integrate such an equality with respect to time. By 
observing that A(^)5t<h + \'(jp)dt^^ = dt{X(jp) $), we simply obtain 

I Qv + V(1 * 0) • Vn + r / 0rnr + « / (1 * ©r) vr = — A(^) <|)n (5.2) 

Jn Jn Jr Jr Jn 

for every {v,vr) E V x Vr such that ur = Hr ™ choosing n = 0 and 

Ur = 0r; integrating with respect to time and forgetting some nonnegative terms on the 
left-hand side, we easily infer that 

||0i ~ 02||L2(o,t;jr) < D ||<hi — for t G [0,T] (5.3) 

where we have marked the constant D for convenience. If we combine 05.31) written for 

= 3^1 (0i) with 05.ip . we deduce the estimate 

||lh(0i) - ^(02)111.(0,= ril(0i - Q2)is)\\j,ds 

Jo 

Jo Jo 

< D‘^C‘^ f II01 - Q 2 \\l 2 ^o,s;H)ds for every t G [0,T] 

Jo 

and this can be iterated. By doing that, we obtain the inequality 

I f LI 

for every 0^ G LJ{Q), t G [0,T] and m > 1. By choosing m such that (C^D^T)'^ < ml 
and t = T, we see that 3”™ is a contraction in L‘^{Q), whence 3” has a unique hxed 
point 0. Then, 0 and the associated functions 0r and $ that enter the construction 
of Ti and T2 provide a unique solution to the linearized problem 02.4ip - 02.44p with the 
regularity 02.38l) - 02.40p . 

Continuous dependence. This is given by the estimate 02.451) we prove at once. By 
integrating 02.411) with respect to time and proceeding as we did for 05.2p . we have 


0n-|- / V(1 *0) ■ Vn-f r / 0rnr + « / (l*0r)'yr = 


/ A(<p)<hn-t-a / m(l*h)vr 

'n Jr 
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for every (p,Pr) E V x Vp such that ur = P|p and a.e. in (0,T). Now, we take v = Q, 
vp = 0r and integrate over (0, t). Besides, we multiply fl2.42p by <59*$, where <5 is a positive 
parameter, and integrate over Qt. Then, we add the equalities we obtain to each other. 
Furthermore, we add the same term {6/{AT)) to both sides for convenience. By 

also owing to the boundedness of the coefficients and to the Young inequality, and setting 
C := ||A(^)|1L, we have 




ler 

f <5 / 


1 

2 




|V(l*0)(f)|' + r 
5(7 




|0rr + 

(5 


a 


^ Qt 


AT 


2 J 
|d>(t)| 


|(l*0r)(t)p 


= — A(<yc)<F0 + a m{l * h) 0i 


^ Qt 


'St 


'Qt 


+ S I [-^{ip) +'d\'{(p))^dt^ + 6 I X{(p)Qdt^ 


<SC / |0|2 + 


'Qt 


A6 


l$r + 


r 


'Qt 

|2 


4T 




|0rr + c 


'St 


'St 


5 

2 


\dtA>\^ + 5c \^\^ + 5C | 0|2 + 


'Qt 


'Qt 


'Qt 


4T 


1 * hp 


Moreover, by arguing as we did to derive the hrst inequality in fl3.18p . we see that 

5 


6 

4T 




\dtA>\ 


^ Qt 


Therefore, by choosing 5 such that 25C < 1 and applying the Gronwall lemma, we obtain 
I|0|U2(Q) + I|l * 0||L°°(O,T;y) + ||0r||L2(s) + ||*h||j^l(o,T;R)nL°°(0,r;U) < c||l * ^I|l2(S) • 


At this point, it is easy to see that the above estimate implies 02.451) . 
Here is the main result of this section. 


□ 


Theorem 5.1. Letu G X. Then, S is Frechet differentiable atu and the Frechet derivative 
[/!§](«) is precisely the map T> G X(X, y) defined in the statement of Proposition 


Proof. We £x n G X and the corresponding state (i?, -dp, ip) := §{u) and, for h G X, we set 

:= S(n + h) and (C^ Cm h'") ^ - 0,4 - 0r, 9?^ - ^ - <h) 

where (0, 0r, *h) is the solution to the linearized problem corresponding to h. We have 
to prove that || (C^) Cm h^)lly/ll^l|x tends to zero as ||h||x tends to zero. More precisely, we 
show that some constant c exists such that 

IKC^) Cm h^)lly ^ '^II^IIl2(s) provided that ||h||x < 1 (5.4) 

and this is even stronger than necessary. So, we assume ||h||oo < 1 and make a preliminary 
observation. As u is hxed and ||m + /i||oo < ||n||oo + 1 for every h under consideration, we 









Colli 


Gilardi 


Marinoschi 


23 


can apply Theorem 12.41 and Corollary 12.61 and find constants eM. and G D{f3) 

independent of h such that 

and •d, < < d* a.e. in Q (5.5) 

^ and <yC. < < </?* a.e. in Q. (5-6) 

Moreover, we can also exploit the second part of Theorem 12.21 with ui = u + h and 
U 2 = u. By doing that and also owing to the continuous embedding V C we derive 

from fl2.19p 

II# - 'd||Loo(o,T;R) < C'||h||i,2(s), \\(p'" - ■^||l4(q) < MC'||h|| l2(s) (5.7) 

IIT’^ ~ <d||L°°(0,r;L4(n)) < MC'WhllL^^T.) (5.8) 

where the constants M and C do not depend on h. Now, the problem solved by (#, #, r]^) 
is the following 


[ dtC^v+ / V#-Vn + r / dtC^vr + a f C^r 
Jn Jn Jr Jr 

= - [ - X{Tp)dtip - X{Tp)dt^ - X'{Tp)dtTp^}v 

Jn 

for every (n,nr) E V x Vr such that vr = n|p and a.e. in (0,T) (5.9) 

dt'T]'^ - aAr]’^ = -E^ + E!^ a.e. in Q, where (5.10) 

El ■■= 7(<d^) - 7(^) - 7'(^) ^ _ (5.11) 

:= #A((^'^) - ^X{Tp) - X{Tp)Q - ^X'{Tp) $ (5.12) 

dnV’^ = 0 a.e. on E and C^(0) = 7^(0) = 0 a.e. in 12. (5.13) 

By integrating fl5.9p with respect to time, we obtain 

[ C’^v+ [ V(1 * #) • Vn + r [ C^vr + a [{1* C^)vr = - [ 

Jn Jn Jr Jr Jn 

for every {v,vr) E V x Vp such that vp = n|p and a.e. in (0,T) (5-14) 

where := ^(#) - 2;(^) - A(^)$ . (5.15) 


At this point, we take n = and ur = Cr 115.14p and integrate over (0,t). Besides, we 
multiply fl5.10p by and integrate over Qt. Finally, we add the resulting equalities to 
each other. We obtain 

[ |V(l*#)#P + r [ |Cr"P + f [ 

Jn JT.t ^ Jr 

+\ I iy(i)t+<^ [ ivyp 

^ Jn JQt 

= - f [ E^7]^+ [ E^r]’^ (5.16) 

and we now estimate each term on the right-hand side. To this end, we represent the 
functions and E’^, i = 1,2, in different forms. By applying the Taylor expansion to A 
and 7, we find functions (p and (f taking values between the ones of Jp and (p^ such that 

W) - Xw) = Mip) [P - i?) + i A'(^) (P - If f 

lip) - iP) = ip) (P + ^ i'P) (i' ~ i)'‘ 
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whence immediately 


In order to rewrite Ag, we observe that, if G : —)■ M is a function, we have for every 

(i/, z) and (?/', z') in and for a suitable (^, z) in between 


G{y\ z') - G{y, z) - \/yzG{y, z) ■ {y -y,z' - z) 


[?/'• 




A DlAyJ)^[y' -y^z' -z]. 


In the particular case G{y,z) = yX{z), the above formula becomes 

y'X{z') - yX{z) - X{z){y' - y) - yX\z){z' - z) 
= \ yX’'{z){z’ - zf + >^{z){y’ - y){z' - z). 


Therefore, there exist functions ■d and Cp taking values between the ones of X} and and 
the ones of Tp and ip^, respectively, such that 

- ^A(^) 

= - d )+- v) + \^\"mp - w + - p 


so that 


^2 = A(^) c" + ^A'(^) + ^ + A'(^)(i9" - - ip). 

Notice that -d, < 'd < 'd* and that p), < p,(p,(p < p*■ Therefore, recalling fl2.35p . the 
right-hand side of fIB.lbp can be estimated by the Young inequality as follows 

_ f _ f ^ f < 1 ^ |^h|2 f |^/a|2 

JQt JQt JQt JQt JQt 

+ c[ / \^^-mp^-p\\y>^\. 

Qt Qt Qt 

<7/ ICl' + c/ \y’^\^ + c [ \p’^-p\^ + c [ \'d’^-mp^-p\\y’^\. 

4 JQt JQt Jq JQt 

On the other hand, thanks to the Holder inequality and to the continuous embedding 
V C T'^(O), we have for every 5 > 0 

[ [)\{^^ -ms)h\\{p^ -ms)\U\\v\s)\Uds 

JQt Jo 

— \\y ('S)llv'^'5 + C 5 llv? ~ ¥^IIl°°(0,T;L'‘(q)) 1119 — 'd||i;,2(o,r;H) • 

Jo 


At this point, we choose 6 small enough, apply the Gronwall lemma and account for 05.711 - 
05.8p in order to estimate the norms of — i? and of p’^ — p. This yields 05.4p and the 
proof is complete. □ 
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6 Necessary optimality conditions 


In this section, we derive the optimality condition fl2.52p stated in Theorem 12.91 We start 
from fl2.33p and first prove 02.461) . 

Proposition 6.1. Letu* he an optimal control and (i9*, i9p, (p*) := S(u*). Then, condition 
02.46P holds. 

Proof. As already said in Section [2], we just have to apply the chain rule for Frechet 
derivatives. Clearly, the Frechet derivative [D'do]{'d,'dT)T) of the functional 3o exists at 
every point of y and it is given by 

[D2iq] (i9, t) ■ (hi, h2, ha) e ^ H- - -dq) hi + K2 [ (^(T) - ifn) h^{T). 

JQ Jn 

In particular, this holds if {'d,'dr,T) = I'd*, hip, (p*) = S{u*). Therefore, Theorem 15.11 and 
the chain rule ensure that 3 is Frechet differentiable at u* and that its Frechet derivative 
[D3]{u*) at any optimal control u* acts as follows 

[D3]{u*) : h e X ^ Ri / (r - i9q) 0 + K 2 / (<^*(T) - p^n) $(T) 

Jq Jn 

where (0, 0r, $) is the solution to the linearized problem corresponding to h. Therefore, 
fl2.46p immediately follows from fl2.33p . □ 


The next step is the proof of Theorem 12.81 For convenience, we consider the equivalent 
forward problem in the unknown {p,pr,(j) given by {p,pr,P)(t) ■= {p,Pr,(l)(T — t) and 
corresponding to the new coefficient and given terms dehned accordingly. However, to 
simplify notations, we write p, pr and q instead of p, pr and q in the sequel. The new 
problem is to find {p,pr,q) satisfying fl2.47p - fl2.48p and solving 


/ dtpv + / Vp ■ Vn + / aqv + t dtPr vr + a prW = 
Jn Jn Jn Jr Jr 

for every v gV, vr := n|p, and a.e. in (0,T) 

/ dtqv + a / Vg ■ Vn + / bqv = / gdtpv 
Jn Jn Jn Jn 

for every n G C and a.e. in (0, T) 

p(0) = 0 and g(0) = go a-e. in 12 


fv 


( 6 . 1 ) 


( 6 . 2 ) 

(6.3) 


where a, /, b, g and go are deduced from fl2.49l) - fl2.51l) . Thus, we have a, b, g E 
f G LP‘{Q) and go G V. We aim to use a contraction argument in a weaker functional 
framework. However, it will be clear from the proof that the unique solution we find 
satishes fl2.47p - fl2.48p . 


The equivalent fixed point problem. For a given g G L‘^{Q), we consider the problem 
obtained by writing flb.ip with g replaced by g and the initial condition p(0) = 0. If 
we introduce the spaces V and TC and the bilinear form a given by fl3.22p - fl3.25p . and 
argue as in the proof of Lemma 13.21 we see that the problem under consideration has a 
unique solution (p,pr) satisfying p G 77^(0, T; 77) D L°°(0, T; V) and fl2.48p . However, p is 
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smoother since one can argue as in Remark 12.11 to derive Ap G L‘^{Q) and dnP G 
from the regularity already achieved. We set 3^i{q) := p and 3^i{q) := {p,Pr)- By doing 
that, we hnd a map Ihi : L‘^{Q) —)■ H) and an associated map 3^i that we use in 

the rest of the proof. Now, for p G H), we consider fl6.2p complemented by the 

second initial condition in fl6.3p . As 6 G L°°{Q), gdtp G L‘^{Q) and go ^ V, such a problem 
has a unique solution q satisfying 02.471) . and we set g := 3^2(p)- We thus obtain a map 
3^2 : 37^(0, T; hr) —)■ L^{Q) and consider 5” : L‘^{Q) —)■ L‘^{Q) dehned by 3^ := 3^2 ° 3^i- 
Let us point out that, for a given g G L'^{Q), 5'(g) actually takes values in the space 
H^{0, T; ih) n L°°(0, T; 1/) fl L^(0, T; W). The problem under consideration is equivalent 
to the existence of a unique hxed point for T. Indeed, if g is such a hxed point, g and 
the corresponding p and pr provide a solution satisfying 02.47^ - 02.48p by construction. 
Conversely, if (p, Pr, g) solves the problem, then g is a hxed point of T. 

The contraction argument. It suffices to hnd a constant C such that 

||(T(gi) -T(g 2 ))(t)||H < C Hg^ - g 2 1^2(0,(6.4) 

for every gi, g 2 G L‘^{Q) and every t G [0,T]. Indeed, this implies that some iterated 
of T is a contraction. In order to prove 06.4p . we take any pair of functions g^ G 
i = 1,2, consider the pairs (pi,Pi,r) •= 3^i(^i) ^md the functions g^ := T(gj) and write 
their dehnitions, i.e., 06.ip written with g^ in place of g and 06.21) with pi in place of p. 
Then we take the two diherences. We set for convenience g := gi — g 2 and similarly 
introduce p, pr and g. At this point, we formally test the hrst diherence by n = dtP 
and fr = dtpr (by avoiding the technical approximation of such test functions by C- and 
l/p-valued functions) and integrate over (0,f). At the same time, we multiply the second 
one by g and integrate over Qt- We have 

f \dtpf + l [ |Vp(t)|^ + r f l^fprl^ + l f |pr(t)r 
jQt ^ Jn JT,t ^ Jv 

as well as 

+ |Vg|^ = - / h\q\^ + f gdtpq. 

J Q. J Qt J Qt ^ Qt 

As a, b and g are bounded functions, it is straightforward to deduce that 

\\p\\m(0,t-,H)nL<=^{0,t-,V) + ||Pr||L2(0,4;Hr) < Cl \\q\\L2{0,f,H) 
ll'?l|L°°(0,t;H)nL2(0,t;y) < C 2 \\p\\h^{0PH) 

for every t G [0,T] and some constants Ci and C 2 independent of g^ and t. By combining 
such inequalities, we deduce that 

\\q\\L°°{Q,t-,H) < C 1 C 2 ||g||L2(0,t;H) 

whence fl6.4p follows with C := C 1 C 2 . Thus, Theorem 12.81 is completely proved. □ 

At this point, we are ready to prove Theorem 12.91 on optimality, i.e., the necessary 
condition fl2.52p for u* to be an optimal control in terms of the solution (p,pr,g) of the 
adjoint problem fl2.49p - fl2.5ip . So, let u* be an optimal control and hx an arbitrary 
u G IXad- We write both the variational formulations of the linearized problem at the 
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optimal state (i9*, i9p, (f*) := S(u*) corresponding to h = u — u* and the adjoint system. 
We have a.e. in (0, T) 


/ dtev+ / V0-Vn+ / <!> v 

Jq Jq Jo Jo 

+ r J dtOr^r + Oi j 0rnr = « Jm{u — u*)vr 

f dt^v + a f V$-Vn+ f ^'{ip*)<!>v= [ rX'{^p*)^v+ f X{ip*)Qv 
Jo Jo Jo Jo Jo 

- / dtpv + / Vp ■ Vn - / X{ip*)qv -r dtpr vr + a prvr 

Jo Jo Jo Jr Jr 

= Ki [ {iX* - dQ)v 

Jo 

- [ dtqv +a f X7q-Vv+ f {y(if*) - iJ*X'(p*))qv = f X{ip*)dtpv 


(6.5) 

( 6 . 6 ) 


(6.7) 

( 6 . 8 ) 


where fl6.5l) and fl6.7l) hold for every (n,nr) G C x Vp such that vr = , while fl6.6p 

and fl6.8l) are required for every v G V. Moreover, the functions at hand satisfy the 
homogeneous initial or hnal conditions as specihed in fl2.44p and fl2.51l) . We choose 
(n,nr) = {p,Pr) in fl6.5p . v = q in fl6.6p . {v,vr) = (—0, —0r) in fl6.7p and v = —$ 
in fl6.8p . Then, we sum all the resulting equalities to each other and integrate over (0,T). 
Several terms cancel out and we obtain 


'Q 


{dtQp + X{ip*)dt^p+ X\ip*)dt(p* <hp + dt^q + dtpQ + dtq^ + X{ip*)dtp^} 

) 

+ T {dtQrPr + dtPr 0r) 


= a m{u — u*)pr — / {J}* —'Jq)Q 


IQ 


As the expression between braces is equal to dt[Qp + A(<p*)<hp + <hg), the above equality 
becomes 


/ (0p + A((p*)$p + $g)(T) - / (0p + A(<p*)$p + $g)(O) 

Jo Jo 

+ T [ (0rPr)(T)-r [ (0rPr)(O) 

As As 

= a m{u — u*)py — ki ('d* — 'Jq)Q ■ 

As Jq 

Owing to the relations 0r(O) = 0(O)|p and Pr(T) = p(T)|p, and accounting for the initial 
conditions fl2.44p and the hnal conditions fl2.5ip . the above equation reduces to 

f ipm n,_{<p‘{T) - <pa) 


= a m{u — u*)pr — Ri / {'&* — 'iJq)Q ■ 

Jt. Jq 
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Therefore, the inequality fl2.46jl we have already established in Proposition 16.11 implies 


a / m(M — u*)pt > 0 for every u G Uad ■ 

is 


Moreover, the positive coefficient a can be removed. At this point, a standard argument 
leads to the pointwise relation fl2.52p and to its consequences listed in the statement, and 
the proof of Theorem 12.91 is complete. 

7 Appendix 

This section is devoted to a rigorous proof of fl3.34p . With respect to the formal procedure 
of Section [21 we replace derivatives with difference quotients, essentially. For h G (0,T) 
and V G L‘^{Q) or n G T^(S), we dehne on (0,T — h) by setting v^{t) := v{t + h). 
We integrate fl3.4p with respect to time over (s, s + h) and test the equality we obtain by 
V = — 'de)(s) and ur = r fhe same time, we write fl3.5p at times s + h 

and s, multiply the difference by — <Pe)(s) and integrate over with respect to space. 
Finally, we add the equalities obtained this way to each other and integrate over (0,f) 
with respect to s. We have 



(7.1) 


All the terms on the left-hand side are nonnegative, the last one by monotonicity, and 
the hrst integral on the right-hand side can be estimated in an obvious way by using the 
Lipschitz continuity of vr. Moreover, we have 



The critical term is the next one, since the cancellation of the formal procedure does not 
occur here. We observe at once that the functions have a common Lipschitz constant 
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since A is Lipschitz continuous. We apply the mean value theorem to Ag and obtain a 
function taking values between the ones of and (p^, in order that 

[ { “ ^eK{^s)) {^e “ {^e “ ^s) } 

jQt 

jQt 

<C [ -l^eW^e -¥’e\'^ + C [ 1^1 

^ Qt ^ Qt 

We treat the last two integrals, separately, by owing to the Holder and Young inequalities 
and to the continuous embedding V C As such embedding is also compact, fl3.30p 

and m Sect. 8, Cor. 4] imply that the functions ip^ are equicontinuous L"^(H)-valued 
functions. Hence, for every 5 > 0, there exists > 0 such that ||</?£(s) — <y5£(s)||4 < 6 
for every e G (0,1) and s G [0,T — h] whenever h < hs- In the sequel, 5 is a positive 
parameter, say 6 G (0,1), whose value is choosen later on, and it is understood that h 
does not exceed the corresponding hs- Therefore, we have 

JQt JO 

<C6 r||(i9,"-i9.)(.)||2||((^,"-^.)(.)lkd5 
Jo 

<6 f\\{^':-^,){s)ryds+c6 f\m-^s){s)rHds 
Jo Jo 

[ \\i^e-^e){s)\\lds + ch^\\dtlJe\\l2^o^T-,H)<^ [ \\{^e - ^e){s)\\l ds + C 

Jo Jo 

where the marked constant C satishes ||n ||4 < C||n||v^ for every n G H and the last 
inequality is justihed by fl3.30l) . Next, we have 

jQt Jo 

<6 f\\{^^,-^,){s)\\lds + Cs fmsWvW^e-^e){s)\\lds 
Jo Jo 

and we observe that the function s i-A (5)11^ is estimated in L^(0, T) by 03.301) . At this 
point, we choose 5 small enough and apply the Gronwall lemma. We obtain 

J Qt ^ ^ ^ St 

+ / \i^s - ^em\^ + f 

J Q, J Qi 

<ch^ + c [ (|V(l*i9£)(h)p +|<^£(h)-(/JoP)+c / |(1 * i9£,r)(h)P (7.2) 

Jn Jr 


for h > 0 small enough and for every t G [0, T — h] and £ G (0,1). Now, we estimate the 
last integrals of 07.2p and show that they have order h^. To this end, we argue rather 
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similarly as we did in deriving fl7.1jl . but we use 'de — 'do and (pe — as test functions. 
Namely, we integrate fl3.4p with respect to time over (0, s), test the equality we obtain 
by T = 'de(s) — 'Oq and vr = 'de,r{.s) — "dolp) and integrate over (0,t) with respect to s. 
Besides, we multiply fl3.5p by and integrate over Qt. Finally, we add the resulting 

equalities to each other and suitably rearrange. We have 


\^e{t) - + 7 : / |V(l*4-l*^9o)(t)P + r / |i9£,r(t)-4|] 


+ 0 / Mt)-^o\^ + a / - (^o)P + / - ^0) 

^ Jn JQt 


' Qt 


= - - 7r(<^o)) - ^ 0 ) 


^ Qt 


^ Qt 


- 'doA£(</?o))(9?£ - ^q) - “ '(^o)} 


+ / m(l*M)(i9,,r-i9o|r)+ / V(1 * 4) • V(4 - 4) 

dst JQt 

+ / {crV<^0 • V(<^e: - + (-/5 £(<^o) - 7r(<^o) -'doAe(<^o))(<^£ - T’o)}- (7.3) 

JQt 


Even though we are interested in taking t = h, it is more convenient to let t vary in 
order to apply some Gronwall-type lemma. So, we assume t G [0, h]. Also in this case, all 
the terms on the left-hand side are nonnegative and the hrst integral on the right-hand 
side can be easily dealt with. Moreover, the next term can be treated as in the above 
argument. Furthermore, we have 


/ m(l -'doip) < / \'de,r-'&o\Tp\ +c \l*u\ 

J J Sf J St 

< f + ch'^\\u\\l= f |i9£,r-i9o|rr + ch^ 

J^t J^t 


Next, we point out that 


/ V(l*i9o)-V(4-i9o) 

JQt 

= / V(1 * do){t) ■ V(1 * - 1 * do){t) - [ Vi9o • V(1 * - 1 * i9o) 

Jfl J Qt 

<] [ |V(l*4-l*i9o)(t)r + /i' / |Vi9or+ / ||Vi9o||H||V(l*i9,-l*i9o)(t)||Hdt. 

4 Jn Jn Jo 

Finally, the last integral of fl7.3p can be written as 


fe{(pe-p^o) where fe :=-aAipo - ^^{ipo) - 7i{ipo) --doX^^ipo) 


' Qt 


whence it is bounded by 


\\fe\\H\\Xe-J^o){t)\\Hdt. 
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On the other hand, we observe that 

ph ph 

/ \\V'&o\\Hdt = h\\V'&o\\H = ch and / \\fs\\H dt = h\\fe\\H < ch 

Jo Jo 

by virtne of 02.91) and 03.121) . Hence, we can apply well-known Gronwall-type ineqnalities. 
Namely we combine, e.g., 0 Lemma A.4, p. 156] and 0 Lemma A.5, p. 157]. By ignoring 
some nonnegative terms on the left-hand side, we conclnde that 

[ (|v(l *4 - 1 <ch^ 

Jn Jr 

for 0 <t < h. In particular, we have 


as well as 


Jn 

< + c [ |V(1 * 'do)(^)P < ch‘^(l + llV'dolll') = 

Jn 


< 2(/(J'||i9e,r(s) < ch^. 

Therefore, the right-hand side of 07.2I) has order hJ and the difference quotients associated 
to the terms of the left-hand side of 07. II) are bounded in the proper norms. Hence, 03.34p 
follows, the bound in L°°(0,T; V) for 'dg being due to = dt{l * Vd^). 
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